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Abstract:- The influence of Hall current, chemical reaction and viscous dissipation on an unsteady MHD flow 
and heat transfer along a porous flat plate with mass transfer is investigated. The transformed governing 
equations in the present study were solved numerically by using Galerkin Finite element method. The effects of 
Hall current, viscous dissipation, chemical reaction over primary velocity u, secondary velocity w, heat transfer 
and mass transfer are analyzed. Effects of magnetic parameter M, Prandtl number Pr, Schmidt number Sc, 
Modified Grashof number Gr for heat transfer and Modified Grashof number Gc for mass transfer on the 
velocity components and temperature is also examined. 
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I. INTRODUCTON 

The combined heat and mass transfer problems with chemical reaction are of importance in many 
processes, and therefore have received considerable amount of attention in recent years. In processes such as 
drying, evaporation at the surface of a water body, energy transfer in a wet cooling tower and flow in a desert 
cooler, the heat and mass transfer occurs simultaneously. Chemical reaction can be codified as either 
homogeneous or heterogeneous processes. A homogeneous reaction is one that occurs uniformly through a 
given phase; in contrast, a heterogeneous reaction takes place in a restricted region or within the boundary of a 
phase. A reaction is said to be the first order reaction if the rate of reaction is directly proportional to the 
concentration itself. In many chemical engineering processes, a chemical reaction between foreign mass and 
fluid occur. Das et.al. [1] Studied the effect of mass transfer on flow past an impulsively started infinite 
vertical plate with constant heat flux and chemical reaction. Muthukumarswamy and Ganesan,P [2,3] discussed, 
flow past an impulsively started infinite plate with constant heat and mass flux and effect of chemical reaction 
on moving isothermal surface with suction. 

In the recent years theoretical study of MHD channel flows has been a subject of great interest due to 
its widespread applications in designed cooling systems with liquid metals, petroleum industry, purification of 
crude oil, polymer technology, centrifugal separation of matter from liquid, MHD generators, pumps, 
accelerators and flow meters. In an ionized gas where the density is low and/or the magnetic field is very strong, 
the conductivity normal to the magnetic field is reduced to the free spiraling of electrons and ions about the 
magnetic lines of force before suffering collision; also, a current is induced in a direction normal to both the 
electric and magnetic fields. The phenomenon, well known in the literature, is called the Hall Effect. The study 
of magnetohydrodynamic flows with Hall currents has important engineering applications in problems of 
magnetohydrodynamic generators and of Hall accelerators as well as in flight magnetohydrodynamics. Sato [4] 
has studied the effect of Hall currents on the steady hydromagnetic flow between two parallel plates. Various 
workers Sherman and Sutton [5], Hossain [6], Ram [7], Pop [8] studied the Hall effects. This effect can be taken 
into account within the range of magnetohydrodynamical approximation. Katagiri [9] studied the steady 
incompressible boundary layer flow past a semi-infinite flat plate in a transverse magnetic field at small 
magnetic Reynolds number considering with the effect of Hall current. Agarwal et.al [10] gave the solution of 
flow and heat transfer of a micro polar fluid over a stretching sheet using finite element technique. Sri Ramulu 
et.al [11] studied the Effect of Hall Current on MHD Flow and Heat Transfer along a Porous Flat Plate with 
Mass Transfer. 

Hossain and Rashid [12] investigated the effect of Hall current on the unsteady free -convection flow of 
a viscous incompressible fluid with mass transfer along a vertical porous plate subjected to a time dependent 
transpiration velocity when the constant magnetic field is applied normal to the flow. And an analysis of heat 
transfer taking dissipation function into account, in boundary layer flow of a hydromagnetic fluid over a 
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stretching sheet in the presence of uniform transverse magnetic field has been given by Lodha and Tak [13], 
Gupta [14] studied the Hall current effects in the steady hydro magnetic flow of an electrically conducting fluid 
past an infinite porous flat plate. Sattar and Maleque [15] studied the effect of viscous dissipation as well as 
joule heating on two-dimensional steady MHD free-convection and mass transfer flow with Hall current of an 
electrically conducting viscous incompressible fluid past an accelerated infinite vertical porous plate with wall 
temperature and concentration. Recently, Bala Siddulu Malga et.al [16] studied the effects of viscous dissipation 
on unsteady free convection and mass transfer boundary layer flow past an accelerated infinite vertical porous 
flat plate with suction. In this paper, we investigated Hall current, viscous dissipation and chemical reaction of 
an unsteady free -convection MHD flow of an electrical conducting viscous incompressible fluid with mass 
transfer along an infinite vertical porous plate. We employed Galerkin finite element method to solve the 
governing partial differential equations. The effects of various important physical parameters on the flow and 
heat, mass transfer are analyzed. 

II. MATHEMATICAL FORMULATION 

An unsteady free-convection flow of an electrically conducting viscous incompressible fluid with mass 
transfer along an infinite vertical porous plate has been considered. The flow is assumed to be in x direction 
which is taken along the plate in upward direction. The y axis is taken to be normal to the direction of the plate. 
At timet > 0 , the temperature and the species concentration at the plate are raised to T w T x ) and C w C x ) 
and thereafter maintained uniform. The level of species concentration is assumed to be very low and hence 
species thermal diffusion and diffusion thermal energy effects can be neglected. A magnetic field of uniform 
strength is assumed to be applied transversely to the porous plate. The magnetic Reynolds number of the flow is 
taken to be small enough so that the induced magnetic field can be neglected. The equation of conservation of 
electric change V.J — 0 gives j y =constant, where J — (j x ,j y , j z ) .It is also assumed that the plate is non- 
conducting. This implies j y — 0 at the plate and hence zero everywhere. When the strength of magnetic field is 
very large the generalized Ohm's law, in the absence of electric field takes the form 
W -^JB = a(n e VB+ — 

Bo en e 

where V is the velocity vector; a, the electric conductivity; fi e , the magnetic permeability; (O e , the electron 

frequency; T e , the electron collision time; e, the electron charge; n e the number density of the electron; and 

P e , the electron pressure. Under the assumption that the electron pressure (for weakly ionized gas), the thermo- 
electric pressure and ion-ship are negligible, equation (1) becomes 

Jx =7— r(mu-w) 
Jz = ( u + mw ) 

>z 1+m 2 v ' 

Where u is the x component of V; w, the z component of V and m (=a) e x e ), the Hall parameter. Within 
this frame work, the equations which govern the flow under the usual Boussinesq approximation are: 

d j - 0 ' ' (3) 

I + V % = V P " («' + mw ) + - r -> + ^ C '~ C '»> < 4 > 

dw' ... dw' d 2 w' aHe B o , • ,-, 

— = v— j , (w -mu') (5) 

dt' dy' dy' L p(l+m 2 ) v J ' 

^+v' — = — — + ^(—) 2 (6) 

df dy' pC p dy' 2 pC p \dy'J K ' 



J + ^JB =a(n e VB+^ r VP e) (1) 



dC , _ d 2 C' 

— + v<>y = D — 

dt dy 

Where g is the acceleration due to gravity; B, the volumetric coefficient of thermal expansion; /?*, the 
coefficient of volume expansion with species concentration; T, the temperature of the fluid within the boundary 
layer; C ', the species concentration; p, fi, v, k, c p are respectively density, viscosity, kinematic viscosity, thermal 
conductivity, specific heat at constant pressure; and D, the chemical molecular diffusivity, Chemical reaction 
parameter^ . The initial and boundary conditions of the problem are 
t < 0: u ' — 0, w ' — 0, T — Too, C = C x for ally \ 

t , >Q (u' = 0,w = 0,T = T w ,C'=C w aty' = 0 | (8) 

[u = 0,w' = 0,7 = T m C'= C x , at y -> oo J 
The non-dimensional quantities introduced in the equations (3) to (7) are 



■ + V dy = D — -*i c (7) 
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t = -A y = y -^,(u,V,w) = ^— — '- 

v v Uq 



(C w —Coo) 



(9) 



Where U 0 is the reference velocity. Therefore, the governing equation in the dimensionless form are 
obtained as 



fly 



du 



du 



- — (u + mw) + G r 0 + G r C 

St Sy Sy 2 1+m 2 v ; r c 



Svv ^ Sw S 2 u> 



St 



Sy Sy 2 



1+m- 



(w - mu) 



se se 

St Sy 



i s 2 e 

P r Sy 2 



■ + 



3£ + y3C I (VC 

St 

Where 



st +V ^-f c J^~ KlC 



(10) 
(11) 
(12) 

(13) 
(14) 



M 



is Magnetic field parameter, P r — is Prandtl number 



p u z 



» 2 o 



is Eckert number 



Sc = - is Schmidt number, Ec = 

D Cp(T w -TJ 

q _ vgP(T w T a ) . s jyj oc jjfj ec j Qj-ashof number for heat transfer 
G c = vgP (C ™ Cc °' ) is Modified Grashof number for mass transfer 



t > 0: 



And the boundary conditions (equation (8)) in the non-dimensional form are 
t < 0: u = 0, w = 0, 7 = 0, C = 0 for all y ' 

u = 0,w = 0, 0 = 1, C = laty = 0 (15) 
u = 0, w = 0, 9 = 0, C = 0 as y -> oo. 
From equation (10), it is seen that v is either constant or a function of time t. Similarly, solutions of equations 
(1 1) to (14) with the boundary conditions equation (15) exist only if 

V = -At~ (16) 
Where X is a non-dimensional transpiration parameter. For suction A > 0 and for blowing X < 0. From equation 
(16), it can be observed that the assumption is valid only for small values of time variable. 



III. METHOD OF SOLUTION 

By applying the Galerkin finite element method for equation (11) over a typical two-noded linear 
element (e) ( y } < y < y k ) is 



( yk N T 

yj 



S 2 u<» SuO' 



yk 



hi 



dy 2 

dN Su (e) 
Sy Sy 



St 



- - i^( ma> + u(e) ) + ^ e + G < c) ] d y 



- N T(_ V ^_^_(JL_\ u ie)+ R *\ 
\ Sy St Vl+m 2 / / 



dy =0 



where R* = (- (^) co + G r 8 + G c c), N = [Nj ,N k ], 0« = g], 

u w = N . 0(e), N . = yjgL iNk= y^i i lie) = yk _ yj = h 



The element equation given by 



J [iV fc iV/ N k N k \t u k\ y [NM N k N k \\Vk\ ay + [N k Nj N k N k \[u k \ ay 
, M \\NjNj NjN k ^ Uj] D JNj], . 

(w)kiv ; . ^JkW-* y k \ dy = 0 



+ 



^ k (S^ +R)dy = 0 



where = 



JV ; . JV ; JV ; JV fc I "ft I T", 1 / m \ W N } N } N k l^, 

Af fc AfJwJ Jv fc JV fc JM Vi+m 2 JLjV fc ^ Af fc AfJLuJ' 



NN NN 



N. 



NN NiNv 



NN, NjN k -\\iL-\ 

J at J J \ andR=R* 
N k Nj N k N k \ [u k \ 



Here the prime and dot denote differentiation with respect to y and t. we obtain 



(17) 
(18) 



(19) 
(20) 
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1(0 L-l lJLuJ 2 L-l lJ Lu^J ^ Vl + m 2 / 6 Ll 2JLwJ 

^ (e)= ^[i zlLl -2-0 



Lu fc 

We write the element equation for the elements y^ < y < y and y ; < y < y i+1 . Assembling these element 
equations, we get 

01 R-i 

Ui + 



l 

7>) 



1 -1 
-1 2 
0 -1 

(—)" 2 

\l+m 2 J 



6 



0 

-1 
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1 
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+ 



-1 
-1 
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6 



1 

0 

-1 



1 

4 
1 



U i+1 

u i+1 



1* 1 



(21) 



Now put row corresponding to the node i to zero, from equation (21) the difference schemes with — h is 



+4Uj + u i+1 ) + 

6 



// M \h V l\ (2 / 2M \h\ 

( v lir^J6-2-ftJ U; - i + U + lrwJ3J U; + 



ff-^r) " + - - r) = R and here R = R* ■ 



Using the Cranck-Nicolson method to the equation (22), we obtain: 
Atu'^l + A 2 u{ +1 + A 3 u[H + /4 4 u{_ 1 + 4 5 uf + 4 6 it/ +1 = R 
Similarly, the equations (12), (13), (14) are becoming as follows: 

+ B 2 w/ +1 + B 3 w{£ + B iW l_ t + B 5 wf + B 6 w{ +1 
C x e\ll + C 2 6( +1 ->- r aJ+i 



B t wL 



+ C-$6{+l + C 4 0/L 1 + C 5 0/ + C 6 8 ] i+1 - R B 
D lC {ll + D 2 cj +1 + D 3 c{£ + D,cU + D 5 c/ + D 6 cj +1 = 0 
The initial and boundary conditions (15) reduce to 

u(i,0) = 0,w(i,0) = O,0(i,O) = 0,C(i,0) = 0 /or aiZ t 
u(0,y) = 0,w(0,)) = 0,0(0,)) = 1,C(0,)) = 1 | 

u(ioo,y) = o,w(ioo,y) = 0,0(100,;) = i,c(ioo,/) = oj /or 71 

^-C--+iG= I )-'+T)A-(i+*+?( i s I )> 
^-G+K^-t-^-H^+'-t-S- 

M-M^)t-;)^ = H(^) + 



(22) 

(23) 

(24) 
(25) 
(26) 

(27) 
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-G 
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\6 3 Vl+m 2 //' 

V 6 Vl+m 2 / 2 6/ 



*-(-*-fc=0"-3-*-(-H;=0+'--T-3 



G-f-T). c « = G+f-T> 

\6 Sc 2 / \6 Sc 2 / 



<i-G--E+T)-ft-e+D- 

1-G-s)-<i-C- + S + T)- 

£>! = (l - 6r ^ + 3pV + K^, D 2 = (4 + 12r j- + 4/^fc), 
D 3 = (l - 6r^- 3pV + /Y 1 fc),D 4 = (l + er^-SpV-^fc), 
D 5 = (4-12r£-4Kifc),D 6 = (l + 6r£ + 3pV - Kifc), 

p _ Mm _ ( U l+1 -Ui-! \ 2 



Here r = 77 where k, h is mesh sizes along y direction and time direction respectively. Index i refers to space 



1+m 2 

k 
h 2 

and j refers to time. The mesh system consists of h=0.1 for velocity profdes and concentration profiles and 
k=0.05 has been considered for computations. In equation (23)-(26), taking i=l(l) n and using initial and 
boundary conditions (15), the following system of equation are obtained. 

AiX^Bi, i = 1,2,3 .... " (28) 

Where -i4 ; 's are matrices of order n and X t and B ; ' s are column matrices having n-components. The 
solution of above system of equations are obtained using Thomas algorithm for velocity, temperature and 
concentration. Also, numerical solutions for these equations are obtained by MATLAB -Program, Galerkin finite 
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element method. The same MATLAB -Program was run with slightly changed values of h and k and no 
significant change was observed in the values of u, w, 6, c. 



To get the physical insight of the problem, the results are discussed through graphs. The numerical 
results are obtained for the velocity profdes of main (primary) flow and cross (secondary) flow and the 
temperature profiles, concentration profiles have been shown graphically for different flow parameters, such as 
magnetic field parameter M, Hall current m, Eckert number Ec, Chemical reaction parameter K± Modified 
Grashof number Gr for heat transfer, Modified Grashof number Gc for mass transfer, Schmidt number Sc, 
transpiration parameter X, and time t. The primary velocity profiles are shown in fig.l (a), fig. 2 (a) and fig.3 (a). 
The Hall current effect is to increases the primary velocity component u, and this is because of the fact that Hall 
Effect, in general reduces the level of Lorentz forces effect on the fluid, whose tendency is to suppress the flow. 
The primary velocity flow is decreases with the increase of magnetic field parameter M and the transpiration 
parameter X. It is also observed that the primary velocity u increases with the time t. It can be seen that the effect 
of Gr, and Gc is to increase the primary velocity profiles u, which is due to with an increase in buoyancy force 
due to temperature differences, the flow is accelerated in the boundary layer. Here the positive of Gr (>0) 
correspond to cooling of the plate by natural convection. The effect of Prandtl number Pr is reduces velocity 
profiles, the reason is that the smaller values of Pr are equivalent to increasing the thermal conductivities, and 
therefore heat is able to diffuse away from the heated surface more rapidly than for higher values of Pr. The 
influence of Schmidt number Sc is to reduce the velocity profiles. It is also noticed that the influence of Eckert 
number Ec, on the velocity profiles is increases the primary velocity profiles u, and designates the ratio of 
kinetic energy of the flow into the boundary layer enthalpy differences. It embodies the conversion of kinetic 
energy to internal energy by work done against the viscous fluid stresses. Hence, greater viscous dissipative heat 
causes a rise in the velocity profiles. 

The secondary velocity profiles w are shown graphically in fig.l (b), fig. 2 (b) and fig. 3 (b). The 
secondary velocity profiles w increase with time t, Hall parameter m and the magnetic field parameter M. But 
decreases owing to increase of transpiration parameter X, it is observed from figures. It also can be observed that 
the secondary velocity profiles w increases as the Grashof number Gr (for heat transfer) and Gc (for mass 
transfer) increases, but it decreases as Prandtl number Pr and Schmidt number Sc increases. It also seen that 
effects of viscous dissipation increases the secondary velocity profiles w. 

From the fig. 5 it is evident that the temperature of the fluid increases with the increase of time t. But it decreases 
with the increase of transpiration parameter X and Prandtl number Pr. The viscous dissipation effect on 
temperature profiles shown in fig. 3(c). The viscous dissipation heat causes a rise in the temperature which is 
shown from the figure. 

Fig. 6 is drawn for the concentration profiles. From the figure it can be seen that the concentrations 
profiles are increases with the increase of time t, while it is decreases with the increase of transpiration 
parameter X and Schmidt number Sc. 

The effect of chemical reaction parameter K\ on primary velocity and secondary velocity profiles, 
temperature and concentration profiles are shown in fig. 4 (a)-(d) respectively. From the fig.4 (a) & (b) it can be 
observed that the effect of K l7 leads to decreases the primary velocity profiles as well as secondary velocity 
profiles. Fig4(c) reviles that an increase of chemical reaction parameter leads to decrease the temperature 
profiles. Fig.4 (d) shows that there is fall in concentration profiles due to the increasing the value of chemical 
reaction parameter K\ is noticed. 



IV. 



RESULTS AND DISCUSSION 
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Fig. 1(a): The Primary velocity profiles 
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Fig.3 (b): The Secondary velocity profdes for different values of Ec 
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Fig.3 (c): Temperature Profdes for different values of Ec 




Fig.4 (a): The Primary Velocity Profde for different values of K x 
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Fig.4(c): The Temperature Profile for different values of 
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